Molecular phylogeny of the family Rhabdiasidae (Nematoda: Rhabditida), with morphology, genetic characterization and mitochondrial genomes of Rhabdias kafunata and R. bufonis.

BACKGROUND: The family Rhabdiasidae (Nematoda: Rhabditida) is a globally distributed group of nematode parasites, with over 110 species parasitic mainly in amphibians and reptiles. However, the systematic position of the family Rhabdiasidae in the order Rhabditida remains unsolved, and the evolutionary relationships among its genera are still unclear. Moreover, the present knowledge of the mitochondrial genomes of rhabdiasids remains limited. METHODS: Two rhabdiasid species: Rhabdias kafunata Sata, Takeuchi & Nakano, 2020 and R. bufonis (Schrank, 1788) collected from the Asiatic toad Bufo gargarizans Cantor (Amphibia: Anura) in China, were identified based on morphology (light and scanning electron microscopy) and molecular characterization (sequencing of the nuclear 28S and ITS regions and mitochondrial cox1 and 12S genes). The complete mitochondrial genomes of R. kafunata and R. bufonis were also sequenced and annotated for the first time. Moreover, phylogenetic analyses based on the amino acid sequences of 12 protein-coding genes (PCGs) of the mitochondrial genomes were performed to clarify the systematic position of the family Rhabdiasidae in the order Rhabditida using maximum likelihood (ML) and Bayesian inference (BI). The phylogenetic analyses based on the 28S + ITS sequences, were also inferred to assess the evolutionary relationships among the genera within Rhabdiasidae. RESULTS: The detailed morphology of the cephalic structures, vulva and eggs in R. kafunata and R. bufonis was revealed using scanning electron microscopy (SEM) for the first time. The characterization of 28S and ITS regions of R. kafunata was reported for the first time. The mitogenomes of R. kafunata and R. bufonis are 15,437 bp and 15,128 bp long, respectively, and both contain 36 genes, including 12 PCGs (missing atp8). Comparative mitogenomics revealed that the gene arrangement of R. kafunata and R. bufonis is different from all of the currently available mitogenomes of nematodes. Phylogenetic analyses based on the ITS + 28S data showed Neoentomelas and Kurilonema as sister lineages, and supported the monophyly of Entomelas, Pneumonema, Serpentirhabdias and Rhabdias. Mitochondrial phylogenomic results supported Rhabdiasidae as a member of the superfamily Rhabditoidea in the suborder Rhabditina, and its occurrance as sister to the family Rhabditidae. CONCLUSIONS: The complete mitochondrial genome of R. kafunata and R. bufonis were reported for the first time, and two new gene arrangements of mitogenomes in Nematoda were revealed. Mitogenomic phylogenetic results indicated that the family Rhabdiasidae is a member of Rhabditoidea in Rhabditina, and is closely related to Rhabditidae. Molecular phylogenies based on the ITS + 28S sequence data supported the validity of Kurilonema, and showed that Kurilonema is sister to Neoentomelas. The present phylogenetic results also indicated that the ancestors of rhabdiasids seem to have initially infected reptiles, then spreading to amphibians.

STEINERNEMA ADAMSI N. SP. (RHABDITIDA: STEINERNEMATIDAE), A NEW ENTOMOPATHOGENIC NEMATODE FROM THAILAND.

A new species of entomopathogenic nematode, Steinernema adamsi n. sp., was recovered from the soil of a longan tree (Dimocarpus sp.) in Mueang Lamphun District, Thailand, using baiting techniques. Upon analysis of the nematode's morphological traits, we found it to be a new species of Steinernema and a member of the Longicaudatum clade. Molecular analyses of the ITS rDNA and D2D3 of 28S rDNA sequences further confirmed that S. adamsi n. sp. is a new species of the Longicaudatum clade, which is closely related to Steinernema guangdongense and Steinernema longicaudam. Using morphometric analysis, the infective juveniles measure between 774.69 and 956.96 µm, males have a size range of 905.44 to 1,281.98 µm, and females are within the range of 1,628.21 to 2,803.64 µm. We also identified the symbiotic bacteria associated with the nematode based on 16S sequences as Xenorhabdus spp. closely related toXenorhabdus griffiniae. Furthermore, we have successfully assessed a cryopreservation method for the long-term preservation of S. adamsi n. sp. Successful cryopreservation of this new species will allow for the longer preservation of its traits and will be valuable for its future use. The discovery of this new species has significant implications for the development of effective biological control agents in Thailand, and our work contributes to our understanding of the diversity and evolution of entomopathogenic nematodes.

Molecular identification of a peroxidase gene controlling body size in the entomopathogenic nematode Steinernema hermaphroditum.

The entomopathogenic nematode Steinernema hermaphroditum was recently rediscovered and is being developed as a genetically tractable experimental system for the study of previously unexplored biology, including parasitism of its insect hosts and mutualism with its bacterial endosymbiont Xenorhabdus griffiniae. Through whole-genome re-sequencing and genetic mapping we have for the first time molecularly identified the gene responsible for a mutationally defined phenotypic locus in an entomopathogenic nematode. In the process we observed an unexpected mutational spectrum following ethyl methansulfonate mutagenesis in this species. We find that the ortholog of the essential Caenorhabditis elegans peroxidase gene skpo-2 controls body size and shape in S. hermaphroditum. We confirmed this identification by generating additional loss-of-function mutations in the gene using CRISPR-Cas9. We propose that the identification of skpo-2 will accelerate gene targeting in other Steinernema entomopathogenic nematodes used commercially in pest control, as skpo-2 is X-linked and males hemizygous for loss of its function can mate, making skpo-2 an easily recognized and maintained marker for use in co-CRISPR.

Kin-recognition and predation shape collective behaviors in the cannibalistic nematode Pristionchus pacificus.

Kin-recognition is observed across diverse species forming an important behavioral adaptation influencing organismal interactions. In many species, the molecular mechanisms involved are difficult to characterize, but in the nematode Pristionchus pacificus molecular components regulating its kin-recognition system have been identified. These determine its predatory behaviors towards other con-specifics which prevents the killing and cannibalization of kin. Importantly, their impact on other interactions including collective behaviors is unknown. Here, we explored a high altitude adapted clade of this species which aggregates abundantly under laboratory conditions, to investigate the influence of the kin-recognition system on their group behaviours. By utilizing pairwise aggregation assays between distinct strains of P. pacificus with differing degrees of genetic relatedness, we observe aggregation between kin but not distantly related strains. In assays between distantly related strains, the aggregation ratio is frequently reduced. Furthermore, abolishing predation behaviors through CRISPR/Cas9 induced mutations in Ppa-nhr-40 result in rival strains successfully aggregating together. Finally, as Caenorhabditis elegans are found naturally occurring with P. pacificus, we also explored aggregation events between these species. Here, aggregates were dominated by P. pacificus with the presence of only a small number of predators proving sufficient to disrupt C. elegans aggregation dynamics. Thus, aggregating strains of P. pacificus preferentially group with kin, revealing competition and nepotism as previously unknown components influencing collective behaviors in nematodes.

Co-option of an Astacin Metalloprotease Is Associated with an Evolutionarily Novel Feeding Morphology in a Predatory Nematode.

Animals consume a wide variety of food sources to adapt to different environments. However, the genetic mechanisms underlying the acquisition of evolutionarily novel feeding morphology remain largely unknown. While the nematode Caenorhabditis elegans feeds on bacteria, the satellite species Pristionchus pacificus exhibits predatory feeding behavior toward other nematodes, which is an evolutionarily novel feeding habit. Here, we found that the astacin metalloprotease Ppa-NAS-6 is required for the predatory killing by P. pacificus. Ppa-nas-6 mutants were defective in predation-associated characteristics, specifically the tooth morphogenesis and tooth movement during predation. Comparison of expression patterns and rescue experiments of nas-6 in P. pacificus and C. elegans suggested that alteration of the spatial expression patterns of NAS-6 may be vital for acquiring predation-related traits. Reporter analysis of the Ppa-nas-6 promoter in C. elegans revealed that the alteration in expression patterns was caused by evolutionary changes in cis- and trans-regulatory elements. This study suggests that the co-option of a metalloprotease is involved in an evolutionarily novel feeding morphology.

Steinernema shori n. sp., a new entomopathogenic nematode (Nematoda: Steinernematidae) from India.

In this study, morphological and molecular features were used to identify a new Steinernema sp. from Chhattisgarh, India. Morphological and molecular features provide evidence for placing the new species into the "bicornutum" clade. The new species is characterized by the following morphological features: infective juveniles with a body length of 587 (494-671) µm; a distance from the anterior end to excretory pore of 46 (43-50) µm; a distance from anterior end to nerve ring of 72 µm (61-85 µm); and E% of 88 (77-97). The first-generation males are characterised by 27 genital papillae and very short spicules, with a length of 61 µm (53-67) µm. The SW% and GS% ratio of S. shori n. sp. are 139 (107-190) and 75 (62-90), respectively. The new species is further characterized by sequences of the internal transcribed spacer and partial 28S regions of the ribosomal DNA. Phylogenetic analyses show that S. shori n. sp. is most closely related to S. abbasi, S. kandii, and S. yirgalemense.

Xenorhabdus aichiensis sp. nov., Xenorhabdus anantnagensis sp. nov., and Xenorhabdus yunnanensis sp. nov., Isolated from Steinernema Entomopathogenic Nematodes.

Three bacterial strains, XENO-2T, XENO-7T, and XENO-10T, isolated from Steinernema entomopathogenic nematodes, were found to represent novel Xenorhabdus species. In this study, we describe these new species by whole-genome and whole-proteome phylogenomic reconstructions, by calculating sequence identity scores using core genome sequences, and by phenotypic characterization. Phylogenomic reconstructions using ribosomal and house-keeping genes, and whole-genome and whole-proteome sequences show that XENO-2T and XENO-10T are closely related to Xenorhabdus japonica DSM 16522T and that XENO-7T is closely related to Xenorhabdus bovienii subsp. africana XENO-1T and to X. bovienii subsp. bovienii T228T. The dDDH values between XENO-2T and XENO-10T and between XENO-2T and X. japonica DSM 16522T are 56.4 and 51.8%, respectively. The dDDH value between XENO-10T and X. japonica DSM 16522T is 53.4%. The dDDH values between XENO-7T and X. bovienii subsp. africana XENO-1T and between XENO-7T and X. bovienii subsp. bovienii T228T are 63.6 and 69.4%, respectively. These dDDH values are below the 70% divergence threshold for prokaryotic species delineation. The newly described species are highly pathogenic to G. mellonella larvae, grow at pH between 5 and 9 (optimum 5-7), at salt concentrations of 1-3% (optimum 1-2%), and temperatures between 20 and 37 °C (optimum 28-30 °C). Biochemical tests such as lysine decarboxylase, ornithine decarboxylase, urease, gelatinase, citrate utilization, indole and acetoin production, and cytochrome oxidase tests allow to differentiate the novel species from their more closely related species. Considering these genetic and phenotypic divergencies, we propose the following new species: Xenorhabdus aichiensis sp. nov. with XENO-7T (= CCM 9233T = CCOS 2024T) as the type strain, Xenorhabdus anantnagensis sp. nov., with XENO-2T (= CCM 9237T = CCOS 2023T) as the type strain, and Xenorhabdus yunnanensis sp. nov., with XENO-10T (= CCM 9322T = CCOS 2071T) as the type strain. Our study contributes to a better understanding of the biodiversity and phylogenetic relationships of entomopathogenic bacteria associated with insect parasitic nematodes.

Host Association and Spatial Proximity Shape but Do Not Constrain Population Structure in the Mutualistic Symbiont Xenorhabdus bovienii.

To what extent are generalist species cohesive evolutionary units rather than a compilation of recently diverged lineages? We examine this question in the context of host specificity and geographic structure in the insect pathogen and nematode mutualist Xenorhabdus bovienii. This bacterial species partners with multiple nematode species across two clades in the genus Steinernema. We sequenced the genomes of 42 X. bovienii strains isolated from four different nematode species and three field sites within a 240-km2 region and compared them to globally available reference genomes. We hypothesized that X. bovienii would comprise several host-specific lineages, such that bacterial and nematode phylogenies would be largely congruent. Alternatively, we hypothesized that spatial proximity might be a dominant signal, as increasing geographic distance might lower shared selective pressures and opportunities for gene flow. We found partial support for both hypotheses. Isolates clustered largely by nematode host species but did not strictly match the nematode phylogeny, indicating that shifts in symbiont associations across nematode species and clades have occurred. Furthermore, both genetic similarity and gene flow decreased with geographic distance across nematode species, suggesting differentiation and constraints on gene flow across both factors, although no absolute barriers to gene flow were observed across the regional isolates. Several genes associated with biotic interactions were found to be undergoing selective sweeps within this regional population. The interactions included several insect toxins and genes implicated in microbial competition. Thus, gene flow maintains cohesiveness across host associations in this symbiont and may facilitate adaptive responses to a multipartite selective environment. IMPORTANCE Microbial populations and species are notoriously hard to delineate. We used a population genomics approach to examine the population structure and the spatial scale of gene flow in Xenorhabdus bovienii, an intriguing species that is both a specialized mutualistic symbiont of nematodes and a broadly virulent insect pathogen. We found a strong signature of nematode host association, as well as evidence for gene flow connecting isolates associated with different nematode host species and collected from distinct study sites. Furthermore, we saw signatures of selective sweeps for genes involved with nematode host associations, insect pathogenicity, and microbial competition. Thus, X. bovienii exemplifies the growing consensus that recombination not only maintains cohesion but can also allow the spread of niche-beneficial alleles.

Xenorhabdus bovienii subsp. africana subsp. nov., isolated from Steinernema africanum entomopathogenic nematodes.

Four Gram-negative bacterial strains isolated from Steinernema africanum entomopathogenic nematodes were biochemically and molecularly characterized to determine their taxonomic position. Results of 16S rRNA gene sequencing indicated that they belong to the class Gammaproteobacteria, family Morganellaceae, genus Xenorhabdus, and that they are conspecific. The average 16S rRNA gene sequence similarity between the newly isolated strains and the type strain of its more closely related species, Xenorhabdus bovienii T228T, is 99.4 %. We therefore selected only one of them, XENO-1T, for further molecular characterization using whole genome-based phylogenetic reconstructions and sequence comparisons. Phylogenetic reconstructions show that XENO-1T is closely related to the type strain of X. bovienii, T228T, and to several other strains that are thought to belong to this species. To clarify their taxonomic identities, we calculated average nucleotide identity (ANI) and digital DNA-DNA hybridization (dDDH) values. We observed that the ANI and dDDH values between XENO-1T and X. bovienii T228T are 96.3 and 71.2 %, respectively, suggesting that XENO-1T represents a novel subspecies within the X. bovienii species. Noteworthy, the dDDH values between XENO-1T and several other X. bovienii strains are between 68.7 and 70.9 % and ANI values are between 95.8 and 96.4 %, which could be interpreted, in some instances, as that XENO-1T represents a new species. Considering that for taxonomic description the genomic sequences of the type strains are compared, and to avoid future taxonomic conflicts, we therefore propose to assign XENO-1T to a new subspecies within X. bovienii. ANI and dDDH values between XENO-1T and any other of the species with validly published names of the genus are lower than 96 and 70 %, respectively, supporting its novel status. Biochemical tests and in silico genomic comparisons show that XENO-1T exhibit a unique physiological profile that differs from all the Xenorhabdus species with validly published names and from their more closely related taxa. Based on this, we propose that strain XENO-1T represents a new subspecies within the X. bovienii species, for which we propose the name X. bovienii subsp. africana subsp. nov, with XENO-1T (=CCM 9244T=CCOS 2015T) as the type strain.

Chitin contributes to the formation of a feeding structure in a predatory nematode.

Some nematode predators and parasites form teeth-like denticles that are histologically different from vertebrate teeth, but their biochemical composition remains elusive. Here, we show a role of chitin in the formation of teeth-like denticles in Pristionchus pacificus, a model system for studying predation and feeding structure plasticity. Pristionchus forms two alternative mouth morphs with one tooth or two teeth, respectively. The P. pacificus genome encodes two chitin synthases, with the highly conserved chs-2 gene being composed of 60 exons forming at least four isoforms. Generating CRISPR-Cas9-based gene knockouts, we found that Ppa-chs-2 mutations that eliminate the chitin-synthase domain are lethal. However, mutations in the C terminus result in viable but teethless worms, with severe malformation of the mouth. Similarly, treatment with the chitin-synthase inhibitor Nikkomycin Z also results in teethless animals. Teethless worms can feed on various bacterial food sources but are incapable of predation. High-resolution transcriptomics revealed that Ppa-chs-2 expression is controlled by the sulfatase-encoding developmental switch Ppa-eud-1. This study indicates a key role of chitin in the formation of teeth-like denticles and the complex feeding apparatus in nematodes.

The modulation effect of the Steinernema carpocapsae - Xenorhabdus nematophila complex on immune-related genes in Drosophila suzukii larvae.

Larvae of the invasive pest Drosophila suzukii are susceptible to the Steinernema carpocapsae - Xenorhabdus nematophila complex and an assessment of the immune-regulatory system activation in this insect was performed to understand the response to the nematode infection. The expressions of 14 immune-related genes of different pathways (Imd, Toll, Jak-STAT, ProPO, JNK, TGF-ß) were analyzed using qRT-PCR to determine variations after nematode penetration (90 min and 4 h) and after bacterial release (14 h). Before the bacteria were present, the nematodes were not recognized by the immune system of the larvae and practically none of the analyzed pathways presented variations when compared with the non-infected larvae. However, after the X. nematophila were released, PGRP-LC was activated leading to the gene upregulation of antimicrobial peptides of both the Toll and Imd pathways. Interestingly, the cellular response was inactive during the infection course as Jak/STAT and pro-phenoloxidase genes remained unresponsive to the presence of both pathogens. These results illustrate how D. suzukii immune pathways responded differently to the nematode and bacteria along the infection course.

Redescription and molecular characterization of Panagrellus ceylonensis (Nematoda, Rhabditida, Panagrolaimidae) from India.

A species of genus Panagrellus was discovered from a wet season form of an oriental common evening brown butterfly Melanitis leda. In this study, a detailed description of Panagrellus ceylonensis is provided including the morphometry, light microscopy and molecular (18S and 28S rDNA genes) studies. Morphological studies on the species agree with original description and characterized by having 12501481m long body in females and 9491305m in males, lateral fields with four longitudinal incisures, lip region continuous and 811 m wide, six offset lips with protruding labial sensilla, neck 124173 m long, excretory pore at the level of basal bulb, vulva post-equatorial (V = 6871), vagina anteriorly orientated with heavily muscled vaginal walls, post-vulval uterine sac 111135 m long or 1.72.6 times as long as the corresponding body diameter, tail conical elongate with an acute terminus in both sexes, spicule 7191 m long, ventrally curved having hooked manubrium and bifurcated lamina tip, lamina ventrally curved with dorsal deflexion at about 60% of spicule length, gubernaculum 2631 m long and well developed. Morphologically, the Indian population of P. ceylonensis does not show a significant difference from the type material of P. ceylonensis in the original description. For molecular studies of this species, the sequence of 18S rDNA is obtained for the first time. Phylogenetic trees based on 18S and 28S rDNA sequences are provided in this study. Additionally, bionomics and global distribution of the species of Panagrellus genus are also discussed. In conclusion, our study provides a comprehensive morphological characterisation and molecular marker sequences of 18S, and 28S genes that can be used to support future taxonomical research on this species and emphasizes the importance of combining molecular data with morphological data to describe the species accurately.

Morphological and molecular profiling of an entomopathogenic nematode Steinernema feltiae: Unlocking its biocontrol potential against vegetable insect pests.

A population of entomopathogenic nematodes, belonging to the Feltiae-clade and labelled J13, was discovered in the agricultural soils of the hilly regions of the Union territory of Jammu and Kashmir, India. Based on morphological, morphometric, and molecular analyses, the nematodes were identified as Steinernema feltiae. The J13 nematode isolate was tested in a laboratory assay for its pathogenicity against six major pests of vegetable crops: Pieris brassicae, Plutella xylostella, Helicoverpa armigera, Agrotis iplison, Trichoplusia ni, and Exelastis atomosa. The morphology of the isolated nematode closely matched the original description, except for the adult females, which had prominent epiptygmata instead of the weakly developed, double-flapped epiptygmata described in the original report. Analysis of the internal transcribed spacer and large subunit rRNA data from the J13 nematodes showed 100% similarity to sequences of the type population, indicating that they are conspecific. The virulence assays revealed that the nematode caused 100% mortality in the tested insect pests within 4872 hours, even at the lowest concentration of 50 infective juveniles per insect. The calculated median lethal concentration varied among the pests, with the lowest number of infective juveniles needed to achieve 50% larval killing being 117 for P. xylostella, 181.74 for P. brassicae, 226.35 for H. armigera, and 202.07 for T. ni at 24 hours post-inoculation. These findings suggest that S. feltiae isolated during the present investigation, may be a viable option for the biocontrol of these insect pests in Kashmir valley, India.

Evolution and Diversity of TGF-ß Pathways are Linked with Novel Developmental and Behavioral Traits.

Transforming growth factor-ß (TGF-ß) signaling is essential for numerous biologic functions. It is a highly conserved pathway found in all metazoans including the nematode Caenorhabditis elegans, which has also been pivotal in identifying many components. Utilizing a comparative evolutionary approach, we explored TGF-ß signaling in nine nematode species and revealed striking variability in TGF-ß gene frequency across the lineage. Of the species analyzed, gene duplications in the DAF-7 pathway appear common with the greatest disparity observed in Pristionchus pacificus. Specifically, multiple paralogues of daf-3, daf-4 and daf-7 were detected. To investigate this additional diversity, we induced mutations in 22 TGF-ß components and generated corresponding double, triple, and quadruple mutants revealing both conservation and diversification in function. Although the DBL-1 pathway regulating body morphology appears highly conserved, the DAF-7 pathway exhibits functional divergence, notably in some aspects of dauer formation. Furthermore, the formation of the phenotypically plastic mouth in P. pacificus is partially influenced through TGF-ß with the strongest effect in Ppa-tag-68. This appears important for numerous processes in P. pacificus but has no known function in C. elegans. Finally, we observe behavioral differences in TGF-ß mutants including in chemosensation and the establishment of the P. pacificus kin-recognition signal. Thus, TGF-ß signaling in nematodes represents a stochastic genetic network capable of generating novel functions through the duplication and deletion of associated genes.

The α2δ Calcium Channel Subunit Accessorily and Independently Affects the Biological Function of Ditylenchus destructor.

The α2δ subunit is a high-voltage activated (HVA) calcium channel (Cav1 and Cav2) auxiliary subunit that increases the density and function of HVA calcium channels in the plasma membrane of mammals. However, its function in plant parasitic nematodes remains unknown. In this study, we cloned the full-length cDNA sequence of the voltage-gated calcium channel (VGCC) α2δ subunit (named DdCavα2δ) in Ditylenchus destructor. We found that DdCavα2δ tends to be expressed in the egg stage, followed by the J3 stage. RNA-DIG in situ hybridization experiments showed that the DdCavα2δ subunit was expressed in the body wall, esophageal gland, uterus, post uterine, and spicules of D. destructor. The in vitro application of RNA interference (RNAi) affected the motility, reproduction, chemotaxis, stylet thrusting, and protein secretion of D. destructor to different degrees by targeting DdCα1D, DdCα1A, and DdCavα2δ in J3 stages, respectively. Based on the results of RNAi experiments, it was hypothesized that L-type VGCC may affect the motility, chemotaxis, and stylet thrusting of D. destructor. Non-L-type VGCC may affect the protein secretion and reproduction of D. destructor. The DdCavα2δ subunit gene also affected the motility, chemotaxis, and reproduction of D. destructor. These findings reveal the independent function of the VGCC α2δ subunit in D. destructor as well as give a theoretical foundation for future research on plant parasitic nematode VGCC.

A new species of the genus Paurodontella Husain & Khan, 1968 (Nematoda: Hexatylina, Sphaerularioidea) with its molecular phylogenetic study.

A new species of the genus Paurodontella, Paurodontella minora n. sp., collected from Alborz Province, Iran, is described and illustrated based on morphological and molecular characters. The new species is characterized by its body length of 393 (350-438) µm and 380 (n = 1) µm in female and male, respectively, 6-7 incisures in lateral field, lip region annulated and continuous with body contour, and total stylet 6.1 (5.5-7.0) µm long. Basal pharyngeal bulb with small posterior extension projecting reaching to the intestine. Excretory pore situated at the level of basal pharyngeal bulb region, no post-uterine sac, conical tail, narrowing to a rounded tip, and rare male with slender tylenchoid spicules and adanal bursa. The new species comes close in morphology and morphometrics to four known species of the genus, namely Paurodontella asymmetrica, Paurodontella balochistanica, Paurodontella densa and Paurodontella niger. In molecular phylogenetic analyses using D2-D3 expansion segments of the large subunit rDNA gene sequence, P. minora n. sp. formed a major clade with species of the genera in the family Sphaerulariidae (Paurodontella, Paurodontoides, Veleshkinema and Sphaerularia) and a sister relation with the members in the families Neotylenchidae and Anguinidae with the same clade support values in Bayesian inference.

Remarks on phylogeny and molecular variations of criconematid species (Nematoda: Criconematidae) with case studies from Vietnam.

The family Criconematidae is a remarkable group of nematodes, containing roughly 600 nominal root-ectoparasitic species, of which many species are known to be significant agricultural pests. Strikingly, our phylogenetic analyses based on 18S, D2-D3 of 28S rRNA, and COI mtDNA sequences of criconematid species, supported by tree topology tests (SH and AU tests), revealed that almost all studied genera, including Criconema, Ogma, Crossonema, Discocriconema, Hemicriconemoides, Criconemoides, Mesocriconema, and Lobocriconema, are not monophyletic groups, a finding that is partly contrary to those of previous studies on these groups. Our results suggest that key morphological characters used in the classification of Criconematidae are the consequence of convergent evolution. It is clear from our studies that the species status of at least 40 sequences of criconematid species from GenBank must be either revised or reconsidered, with analyses based on a polyphasic approach that includes different tree- and distance-based molecular species-delimitation methods (bPTP, GMYC, ABGD1, and ABGD2). Our studies found the ABGD2 output of the automatic barcode method to agree remarkably well with established species delimitations, while in general, the four species-delimitation results corresponding to three barcode regions forwarded significantly more putative species compared to those originally considered. This study also characterised for the first time the populations of Criconemoides myungsugae and Discocriconemella hensungica associated with Vietnamese ginseng, one of the most precious and rare ginseng varieties in the world. Although these populations are morphologically in agreement with the original descriptions of C. myungsugae and D. hengsungica, their molecular data display notable variations compared to the sequences deposited in GenBank. These species demonstrate clearly the immense molecular variations that can be observed in several species of the family Criconematidae.

The improved genome of the nematode Parapristionchus giblindavisi provides insights into lineage-specific gene family evolution.

Nematodes such as Caenorhabditis elegans and Pristionchus pacificus are extremely successful model organisms for comparative biology. Several studies have shown that phenotypic novelty but also conserved processes are controlled by taxon-restricted genes. To trace back the evolution of such new or rapidly evolving genes, a robust phylogenomic framework is indispensable. Here, we present an improved version of the genome of Parapristionchus giblindavisi which is the only known member of the sister group of Pristionchus. Relative to the previous short-read assembly, the new genome is based on long reads and displays higher levels of contiguity, completeness, and correctness. Specifically, the number of contigs dropped from over 7,303 to 735 resulting in an N50 increase from 112 to 791 kb. We made use of the new genome to revisit the evolution of multiple gene families. This revealed Pristionchus-specific expansions of several environmentally responsive gene families and a Pristionchus-specific loss of the de novo purine biosynthesis pathway. Focusing on the evolution of sulfatases and sulfotransferases, which control the mouth form plasticity in P. pacificus, reveals differences in copy number and genomic configurations between the genera Pristionchus and Parapristionchus. Altogether, this demonstrates the utility of the P. giblindavisi genome to date and polarizes lineage-specific patterns.

Analysis of repeat elements in the Pristionchus pacificus genome reveals an ancient invasion by horizontally transferred transposons.

BACKGROUND: Repetitive sequences and mobile elements make up considerable fractions of individual genomes. While transposition events can be detrimental for organismal fitness, repetitive sequences form an enormous reservoir for molecular innovation. In this study, we aim to add repetitive elements to the annotation of the Pristionchus pacificus genome and assess their impact on novel gene formation. RESULTS: Different computational approaches define up to 24% of the P. pacificus genome as repetitive sequences. While retroelements are more frequently found at the chromosome arms, DNA transposons are distributed more evenly. We found multiple DNA transposons, as well as LTR and LINE elements with abundant evidence of expression as single-exon transcripts. When testing whether transposons disproportionately contribute towards new gene formation, we found that roughly 10-20% of genes across all age classes overlap transposable elements with the strongest trend being an enrichment of low complexity regions among the oldest genes. Finally, we characterized a horizontal gene transfer of Zisupton elements into diplogastrid nematodes. These DNA transposons invaded nematodes from eukaryotic donor species and experienced a recent burst of activity in the P. pacificus lineage. CONCLUSIONS: The comprehensive annotation of repetitive elements in the P. pacificus genome builds a resource for future functional genomic analyses as well as for more detailed investigations of molecular innovations.

Sequencing and phylogenomics of the complete mitochondrial genome of Allodiplogaster sp. (Rhabditida: Diplogasteridae): A new gene order and its phylogenetic implications.

Gene order has been utilized as a phylogenetic signal for many taxa. However, its phylogenetic performance has not been evaluated in Nematoda. As there is only one nematode mitogenome available to date, in the Diplogasteridae family, we sequenced the mitogenome of Allodiplogaster sp. and constructed a phylogeny for Nematoda using this updated mitogenome dataset. We then compared this phylogeny to one constructed using gene order. The complete mitochondrial genome of Allodiplogaster sp. was 13,953 bp in size and included 22 tRNAs, two rRNAs, and 12 protein-coding genes. To assess how Allodiplogaster sp. is related to other nematode species, we used Bayesian inference and maximum likelihood algorithms to construct phylogenetic trees of the Nematoda. We found that: 1) The target species Allodiplogaster sp. is closely related to Allodiplogaster sudhausi. The topology of the mitogenome based phylogeny was nearly identical to previous phylogenies created using 18S rRNA data, except for the placement of the Strongyloididae family. 2) The maximum likelihood tree constructed using gene order was roughly consistent with the mitogenome-based tree at the family level, but not at the species level. 3) Protein-coding genes were ordered differently in Allodiplogaster sp. versus Allodiplogaster sudhausi; this represents the first report of such a reordering in the class Chromadorea in our study. Our study confirms that gene order represents useful phylogenetic information for the Nematoda: the maximum likelihood tree based on gene order provided additional support for the nematode phylogeny constructed using molecular data.